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Abstract 

We investigate the nucleon's transverse spin-dependent structure function g2(x,Q 2 ) in 
the framework of the operator product expansion and the renormalization group. We 
construct the complete set of the twist-3 operators for the flavor singlet channel, and 
give the relations among them. We develop an efficient, covariant approach to derive the 
anomalous dimension matrix of the twist-3 singlet operators by computing the off-shell 
Green's functions. As an application, we investigate the Q 2 -evolution of g2(x,Q 2 ) for the 
lowest moment case, and discuss its experimental implication. 

1 Introduction 

The nucleon spin structure observed in deep inelastic scattering (DIS) is described by the two 
structure functions g±(x, Q 2 ) and g2(x,Q 2 ). g\ and g 2 are present for the longitudinally and 
transversely polarized targets, respectively. In the framework of the operator product expansion, 
the leading contribution for g\ comes only from the twist-2 operators, while for g 2 the twist-3 
operators also contribute in the leading order of 1/Q 2 in addition to the twist-2 operators [p], 
0]. g2(x,Q 2 ) is unique, and plays a distinguished role in spin physicsf!], ||, §]. Firstly, it is 
a "measurable" higher-twist structure function: In general, it is difficult to extract the higher- 
twist structure functions by experiments because they usually constitute small corrections to the 
leading twist-2 term. However, g2(x,Q 2 ) is somewhat immune to this difficulty, and contributes 
as a leading term to the asymmetry in the DIS using the transversely polarized target. Secondly, 
g2(x, Q 2 ) is related to an interesting sum rule, referred to as the Burkhardt-Cottingham sum rule: 
fodxg 2 (x,Q 2 ) = 0. Thirdly, g 2 {x, Q 2 ) contains information inaccessible by the more familiar spin 
structure function gi(x,Q 2 ): g2{x,Q 2 ) is related to the nucleon's transverse polarization and to 
the twist-3 operators describing the quark-gluon correlation in the nucleon. 

Recently, the first data of g2(x,Q 2 ) have been reported||. Extensive study is expected to 
be performed at DESY, CERN, and SLAC. These future measurements of g2{x,Q 2 ) will give 
the first opportunity to obtain new information about the dynamics of QCD and the nucleon 
structure beyond those from the twist-2 structure functions. Theoretically, the determination 
of the Q 2 -dependence of g2(x,Q 2 ) is very important to extract these physical information from 
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experimental data. Furthermore, the comparison of the Q 2 -evolution itself between theory and 
experiment will provide a deeper test of QCD beyond the conventional twist-2 level. 

The (Revolution of g2{x,Q 2 ) is governed by the anomalous dimensions, which are deter- 
mined by the renormalization of the relevant twist-3 operators as well as the twist-2 ones. A 
characteristic feature of the higher-twist operators is the occurrence of the complicated operator 
mixing under renormalization: Many gauge-invariant operators, the number of which increases 
with spin (moment of the structure function), mix with each other. [] Furthermore, the operators 
which are proportional to the equation of motion ("EOM operators"), as well as the ones which 
are not gauge invariant ("alien operators"), also participate through the renormalization mixing 
11- 

There have been a lot of works on the (Revolution of g2{x,Q 2 )- Most of them discussed 
the flavor nonsinglet case |7|, [10, DJ. Only a few works treated the singlet case||: Bukhvostov, 



Kuraev and Lipatov derived evolution equations for the twist-3 quasi-partonic operators. Re- 
cently, Miiller computed evolution kernel based on the nonlocal light-ray operator technique, and 
obtained the identical results. However, both of these two works employ a similar framework 
based on the renormalization of the nonlocal operators in the (light-like) axial gauge. Balitsky 
and Braun also treated the nonlocal operators although they employed the background field 
method. On the other hand, covariant approach based on the local composite operators is miss- 
ing. Furthermore, some subtle infrared problem occurring in the renormalization of the generic 
flavor singlet operators has been emphasized in |J . Therefore, the computation of the anomalous 
dimensions for the flavor singlet part in a covariant and fully consistent scheme is desirable and 
should provide useful framework. 

In this work, we develop a covariant framework to investigate the flavor singlet part of 
92{x, Q 2 ) based on the operator product expansion (OPE) and the renormalization group (RG), 



by extending our recent work on the flavor nonsinglet part|K| |TTJ. Sect. 2 is devoted to a 
detailed OPE analysis for the flavor singlet part of g2{x,Q 2 ). We list up all relevant twist-3 
flavor singlet operators appearing in QCD. We derive the relations satisfied by these operators. 
In particular, we obtain a new operator identity, which relates the gluon bilinear operator with 
the trilinear ones. Based on these developments, we give a basis of the independent operators for 
the renormalization. In sect. 3, we discuss a general framework to perform the renormalization 
of the twist-3 flavor singlet operators in a covariant gauge. We compute the off-shell Green's 
functions inserting the relevant local composite operators. Infrared cut-off is provided by the 
external off-shell momenta. In this case, the EOM operators as well as the alien operators should 



be included as independent operators. We introduce convenient projection technique [|I2|, |TT| to 
avoid the complexity stemmed from the renormalization mixing of many gauge-noninvariant 
operators. Sect. 4 contains an application of our framework to investigate the (Revolution of 
the lowest (n = 3) moment. We discuss experimental as well as theoretical implication of our 
results. The final sect. 5 concludes with a brief summary. 



2 OPE analysis of twist-3 flavor singlet operators for g<i{%, Q 2 ) 

Let us start from a quick look at the factorization theorem in QCD[^, [13]. As is well known, 
the cross section for the DIS is given by the "cut diagram" corresponding to the discontinuity 
of the forward virtual Compton amplitude between the virtual photon with the momentum 

*Recently, it has been proved that the twist-3 structure functions obey simple DGLAP evolution equation for 
-/V c — > oo, and that this simplification is universal for all twist-3 nonsinglet structure functions Q. 
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(q 2 = —Q 2 ) and the nucleoli with the momentum P^ (P 2 = M 2 with M the nucleon mass). 
In general kinematics, this virtual Compton amplitude contains all the complicated interactions 
between the virtual photon and the nucleon, possibly including the "soft interactions" where the 
soft momenta are exchanged. However, drastic simplification occurs if one goes to the Bjorken 
limit Q 2 — > oo with x = Q 2 /2P ■ q fixed: The amplitude is dominated by the contribution 
which is factorized into the hard (short distance) and the soft (long distance) parts, and the 
other complicated contributions are suppressed by the powers of l/Q. The factorized amplitude 
corresponds to the process where a parton carrying the momentum k = £P (0 < £ < 1) comes 
out from the soft part, followed by the hard hitting by the virtual photon, and then goes back to 
the soft part. For the case of the polarized target, the structure functions gi(x, Q 2 ) and g 2 (x, Q 2 ) 
appear in the cross section, and they are given by, corresponding to the factorized amplitude, 

gi (x,Q 2 ) = Y,[jm^ 2 )Hu(^Xa s {^, (1) 

g 2 (x,Q 2 ) = E£f^ 2 (e^ 2 )^(|^,« s (^. (2) 

where <f)\,4> 2 are the parton distribution functions corresponding to the soft part. As is well 
known, the twist- 2 distribution 4>\(£, is interpreted as the probability density to find a parton 
of type i (= gluon, u, % d,d,- ■ •) in the nucleon, carrying a fraction £ of the nucleon's momentum. 
The summation of ([]]), (fj) is over all the possible types of parton, i. Hu,H 2i denote the hard 
parts; they correspond to the hard scattering coefficients for the scattering between the virtual 
photon and a parton i. They are calculable systematically by perturbation theory, and the 
dependence on the strong coupling constant a s = g 2 /^n is explicitly shown. In ([!]), (0), the soft 
and the hard parts are divided at the renormalization scale \x. 

The soft part (parton distribution function) is a process-independent quantity which appears 



universally in various hard processes. Following [rj|] this is written as Fourier transform of 
nucleon matrix elements of nonlocal light-cone operators, reflecting the light-cone dominance in 
hard processes. For our case, the relevant quantity is 

$;(^ 2 ) = P + j ^e^{PS\^MlM^AUz)\PS). (3) 

Here z is a light-like vector z 2 = 0, \PS) is the nucleon state with momentum P and spin S 
(P ■ S = 0, S 2 = —M 2 ), and ipi is the quark field. Our Lorentz frame is chosen as P ■ z = P + z~ . 
The bilocal operator is renormalized at the scale \i. [0, z] — Pexp ig ctz) is the 

path-ordered gauge factor, and therefore (H) is gauge- invariant. We considered the case where 
the quark of flavor i comes out of the nucleon, followed by the coupling with the vertex 7^75, and 
then goes back to the nucleon. The quark has £P + as the plus component of the momentum. 

In principle, one can insert various gamma matrices in between the quark fields. In fact, 
one can generate all quark distribution functions up to twist-4 by inserting all the possible 
gamma matrices [0]. For example, if 7 M were inserted in place of 7^75, we would obtain the 
spin-independent distribution functions of twist-2 and twist-4, corresponding to the structure 
functions Fi(x,Q 2 ), F 2 (x,Q 2 ). In our case of (§), additional 75 distinguishes the helicity of 
the quark. Therefore, ([$]) gives the spin-dependent distribution functions <j)\, <ft 2 relevant to the 



structure functions g\(x, Q 2 ) and g 2 (x,Q 2 ). Explicitly, we havefll 



(4) 
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up to the twist-4 corrections. Here Sj_ M = — P^S ■ z)/(P ■ z) + M 2 z^(S ■ z)/(P ■ z) 2 is the 
projection of the spin vector into transverse direction. By comparing both sides of (01), one can 
obtain explicit operator definitions of the relevant distribution functions. 

For further development, it is convenient to employ the moment given by J d££ n -1 3>* (£, /i 2 ), 
which is related to the moments of the structure functions M n [g 1 (Q 2 )} = Jq dxx n ~ 1 gi(x,Q 2 ), 
M n [g 2 (Q 2 )] = Jq dxx n ~ 1 g 2 (x,Q 2 ) (see ([!]), (fj)): As seen from (|3|), going over to the moment 
space corresponds to expanding the nonlocal operator at small quark- ant iquark separation. This 
gives the OPE. We obtain (we consider the case of odd n which is relevant for the DIS) 

J\e-%(t^) = jjVr -1 /i 2 ) + <(^ 2 

= A m • " " A^(PS\^(0) W D^ ■ .-iD^MtylPS), (5) 

where = j[P ■ z), and = <9 M — igA^ is the covariant derivative. The local composite 
operators in the r.h.s. are symmetric and traceless for the indices fii, ■ • ■ , /i n _i, but have one free 
index fi. If one recalls that the twist of the local composite operator is defined as "dimension 
minus spin" , one immediately concludes the following: The symmetrization of the index fi with 
the other indices and the subtraction of the trace terms corresponding to g mh give the twist-2 
operators; these twist-2 operators correspond to gi(x,Q 2 ). On the other hand, the antisym- 
metrization of /x with the other indices introduces one unit of twist, giving the twist-3 operators; 
the matrix element of these operators, having antisymmetric pair of indices, should involve S^^, 
and therefore contributes to g 2 (x,Q 2 ) (see (f|)). 

From (U) and (|j), one sees that the twist-2 operators also contribute to g 2 (x, Q 2 ). As is well 
known, these contributions can be conveniently extracted by@ 

92(x,Q 2 ) = - gi (x,Q 2 )+ f l dy 9l{v ' Q2) +g 2 (*,Q 2 ), (6) 

J x y 

where we denote the twist-3 part by g 2 (x, Q 2 ). Now M n [g 2 (Q 2 )] is given by matrix element of 

2(n — 1 "l r_ i 

^w-Mn-! = $n+ijs LsAS UYi^ 1 ■ ■ ■ D^tfj - traces, (7) 

n 1 J 

where S symmetrizes /ii, • ■ • , /i„_i, and A antisymmetrizes a and /ii as Af"^ 1 = Mf 7 ^ 1 — f^ ia )- 
"— traces" stands for the subtraction of the trace terms to make the operators traceless. n = 
3, 5, 7, ■ • -, because only the charge-conjugation even operators contribute to the DIS. Here and 
below we suppress the flavor index of the quark field: For the flavor nonsinglet part, the flavor 
matrices Aj should be inserted between the quark fields. On the other hand, for the present case 
of the singlet part, the flavor indices have to be summed over. 

In addition to (0), one can construct other types of gauge-invariant operators of twist-3 || [?]]: 

RZ 1 '"^- 1 = ^-(Vi-V n - 1 - l + U l + U n - 1 -i) (i = l,...,n-2) , (8) 
2n 

R^-^- 1 = i n Smi)Y^D^ ■ ■ ■ D^j^- 1 ^ - traces , (9) 



r? <T^i--^„_i -n™ l c 

R n ,E = 1 — — <-> 



where 



2n 

+ ^7 <T 7 5 D /X1 ■ • ■ D^- 2 ^' 1 (ifl- m)ip] - traces , (10) 

Vi = i n gS^ b D^ ---G^ 1 ■■■D^- 2 ^ n - 1 ^ -traces, (11) 

Ui = i n ~ 3 gS^D fll ---G' 7f11 ■■■D fln - 2 Y n - 1 ^- traces. (12) 
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Here m represents the quark mass (matrix). The combination in the r.h.s. of (§) constitutes an 
independent set of charge-conjugation even operators. The operators in (B contain the gluon 
field strength G^ y or the dual tensor G^ u = ^s^a/sG" 13 explicitly; this implies that they represent 
the effect of quark-gluon correlations. Rn,E is the EOM operator; it vanishes by the use of the 
QCD equations of motion, although it is in general a nonzero operator due to quantum effects. 
Rn, m is due to the quark mass effect. 

The above operators are not all independent. By using = {•y fJi ,ft}/2 and [D^D,,] = 
—igG^, it is straightforward to derive the following equation || 

n-1 n ~ 2 

RZi'^- 1 = Kfm" M " -1 + £(n - 1 - W" Mn_1 + KT^ 1 ■ (13) 

n 1=1 

Therefore we can exclude one operator among (|7|)-(|T0|) to form an independent basis. A conve- 
nient choice of the independent operators will be (||), ® and (|T0l). 

For the flavor nonsinglet part of g2{x,Q 2 ), the operators (|8])-(|T0|) form the complete set of 
the twist-3 gauge- invariant operators. For the singlet part, however, this is not the whole story: 
The QCD radiative corrections could replace the quark fields involved in (|7|), (||), ( [TT]) and ([12| ) 
by the gluon fields, and generate the operators which are bilinear or trilinear in the gluon fields. 
In addition, because the QCD lagrangian in a covariant gauge involves the ghost field as a 
dynamical variable, the alien operators involving these unphysical degrees of freedom could also 
participate. These considerations lead to the following "new" operators: 



1 n,G 






A S [G U ^D 


^2 . . . JJ^n- 


X G° —traces, 






(14) 


rpCTfMl- 

1 n,l 




= 






■ ■ ■ G^ H ■ ■ 






- traces (I 


= 2, 


n — 1 . , 


1 n,l 


••/ln-1 




2 gS 






■ D^Gf- 1 ' 




- traces (I 


= 2, 


n — 1 . , i 
■■-.—). ( 16 ) 


1 n,B 




= 


l S{G CT ^ 1 D^ 2 ■ ■ ■ 




—d^-HdTAl) 
a 


+ gf abc {d^ 




)X C } - traces(17) 


rpCTfll- 

1 n,E 


"Mn-1 


= 


-iS{G*Pi D P2 ... 


D^- 2 } a {(D u GJ Mn - 1 ) a 











+gipt a ^- 1 ij + -d^- 1 (d v A a v ) - gf abc (d^-^)x c ] - traces. (18) 
a 



Here the gluon field A^ and the covariant derivative are in the adjoint representation, x and 
X are the ghost fields, and a is the gauge parameter. t a is the color matrix as [t a ,t b ] = if abc t c , 
Tr(t a t b ) = \5 ab . T n> i and T n> i are trilinear in the gluon field strength and its dual tensor, and thus 
represent the effect of three gluon correlations. In contrast to the quark-gluon operators (§), (|H|) 
and flT2]), T n> i (T n> i) is charge-conjugation even by itself. T n ^E is the EOM operator, and vanishes 
by the naive use of the equations of motion for the gluon. T n B is the BRST-exact operator j5] 
which is the BRST variation of the operator: z n_1 «S{G°'^ 1 I? w • • • D^- 2 } 01 ^- 1 ^ - traces. 

Again, the operators (|14])-(|T8"1) are not all independent: Firstly, due to the Bose statistics 
of the gluons, we obtain the symmetry relations T^ 1 = T^Z' l IJ ' n ~ 1 and T^ 1 Mn_1 = 

"^nn-r™" 1 ' Therefore, we can choose T njZ and T n i for I = 2, • • • , as independent oper- 
ators, as indicated in fll5|) and (]T6|) . Secondly, one can derive the relation between T n ^ and 
T n ,i- 

j^O-tll-'-Hn-l ST^ (QU—2—j G^~ 2 \ ( \\ l ' r p a ^ 1 ■■■fJ-n-l _|_ lQTl—2—i C ]\ jrpV^l-- 

n,j / j y / 2 / 2 J \ J n,l y j 2 J v / n,j 

1=2 

5 



(n-l)/2 



+ E ct" 



~in—2—j 



c n n zlzi) (-i)XT 



IrpcrfJtf-fhi-i 



(19) 



where = n!/[r!(n — r)!]. This relation is also due to the Bose statistics of the gluons. Thirdly, 
the analogue of (|Hf), relating the gluon bilinear operator with the trilinear ones, can be obtained 

as 



rp^Hl—Hn—l 
1 n,G 



n-1 

2 

E 

1=2 L 



(1-2) 



c n-2 _ 



2n 



XI n,l 



(n-l-2) 



2 n 



-a 



n-2 
n—l 



n(n 



21) 



2(n-2)(n- 1) 



n-2 



n-2 
l-l 



:-i) 



i+i 



n 



rpCrfJ,l--U n -l 

' 1 n,E 



n 



Note that the quark-gluon-quark operator R n> i of (^]) as well as the operators flTSp, ( |TTD and 
appears in the r.h.s. To derive this relation, we have used [D^, D v ] = —igG^ v and the identities: 

D a G va -\- D p G aa -\- D a G ' av 0, D a G va -\- D v G acr -\- D a G av s vaa pD \G ^, (21) 

where the first identity is the usual Bianchi identity while the second one is a consequence of 
the relation g^e^s = g m e v p^ 5 + g^au-yS + 9 m e a p vS + g^ap-yv The operator identity © is 
new, and is one of the main results of this work. As a result of (|19|), (PI), we can conveniently 
choose a set of independent operators as (|15D , (|17D and flT8|) . 

To summarize our independent basis of the twist-3 flavor singlet operators for the n-th mo- 
ment (see (||)-(0), flTop, (|T7|) and (0)): The quark-gluon-quark operators R n j (/ = 1, • • • , n — 2); 
the quark EOM operator R n> E] the quark mass operator R n<m ] the three-gluon operators T n> i 
{I = 2, • • • , ^y-); the gluon EOM operator T n ^\ the BRST-exact operator T ni #. Among them, 
R n ,i, R n ,E, R n ,m, and T Ui i are gauge- invariant, while T Ut E and T n> B are not gauge-invariant but 
BRST-invariant. These (3n + l)/2 operators will mix with each other under renormalization. 



3 Renormalization of twist-3 singlet operators 

The Q 2 -evolution of the flavor singlet part of g2{x,Q 2 ) is governed by the anomalous dimen- 
sions, which enter into the RG equation for the relevant twist-3 singlet operators obtained in 
sect. 2. In this section we discuss the renormalization of these operators to obtain the anomalous 
dimensions. 

We follow the standard method to renormalize the local composite operators ||. We multiply 
the operators discussed in sect. 2 by a light-like vector A Mi to symmetrize the Lorentz indices 
and to eliminate the trace terms: A w • • ■ A ftn _ 1 E^'" t * n - 1 = A • R a np A Ml • • • A^T^ 1 "^" 1 = 
A ■ T° h etc. We then embed the operators Oj = A • R° t , A • T° u etc. into the three-point 
function as (Q\TOj(tt)A^{x)^(y)ii(z)\ti), (0|TO i (0)A AI (x)X(j/)A p (/)|0), etc., and compute the 
1-loop corrections. We employ the Feynman gauge (a = 1) and renormalize the operators 
in the MS scheme. To perform the renormalization in a consistent manner without subtle 
infrared singularities Q, we keep the quark and gluon external lines off-shell; in this case the 
EOM operators as well as the BRST-exact operators mix through renormalization as nonzero 
operators. 

One serious problem in the calculation is the mixing of the many gauge-noninvariant as well 



as BRST-noninvariant EOM operators. As explained in |ID|, 115 1, gauge-noninvariant EOM op- 



erators are given by replacing some of the covariant derivatives by the ordinary derivatives 
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<9 Mi in ([Top . Furthermore, in the present case, the BRST-noninvariant EOM operators obtained 
similarly from the gluon EOM operator (113) will also participate in the mixing. However, the 



problem can be overcome by direct generalization of the method employed in [11 , |12|| . We intro- 
duce the vector Q r (r = 1,2,3) satisfying A^f2^ = for each external gluon line, and contract 
the Green's functions as fiJ(0|TO i (0)A'*(ar)^(2/)^(2)|0), nlnln 3 p (0\TO j (0)A^(x)A"(y)AP(z)\0), 
etc. This brings the two merits: Firstly, the tree vertices of the gauge (BRST) invariant and the 
gauge (BRST) noninvariant EOM operators coincide. Thus, essentially only one quark (gluon) 
EOM operator is now involved in the operator mixing. Secondly, the structure of the vertices 
for the twist-3 operators are simplified extremely, and the computation becomes more tractable. 

We perform the 1-loop calculation for the one-particle-irreducible three-point function with 
the insertion of A • R° t , A • T° h etc. Fig.|l] shows the relevant 1-loop diagrams. 



A 

(d) 





(i) 0) (k) 



Figure 1: One- loop diagrams (full lines: quarks; wavy lines: gluons) 
The results can be summarized in the following matrix form: 



/ A 


per 


\ 




fz[ k F 


yFF 

im 


yFG 


yFG 


yFF 
L iE 


yFG 
^iE 




( A 






A 


per 

- rt n,m 









yFF 

mm 
















A 


ptT 




A 


rpa 
1 n,l 






yGF 
A lk 


yGF 


yGG 


yGG 


yGF 
^lE 


yGG 
^lE 




A 


rpa 

n,j 




A 


rpa 

1 n,B 















yGG 
^BB 


yGF 
^BE 


yGG 
^BE 




A 


rpa 
1 n,B 




A 




















yFF 
^EE 


yFG 
^EE 




A 


Rn,E 




V A 


rpa 

1 n,E 


) 


bare 


V o 











yGF 
^EE 


yGG 
^EE 


) 


V A 


1 n,E 


) 



(22) 



where the operators with (without) the suffix "6are" are the bare (renormalized) quantities. In 
the MS scheme we express the renormalization constants Z~fi as 



7 AA< 

j pp' 



= 5 pp i5aa' + 



9 



-v-AA' 



(23) 



n 



with D the space-time dimension. Here p = 1, • • • , n — 2, m, E for A = F, while p = 2, 
l)/2, B,E for A = G. The renormalization constant matrix is triangular because the physical 
matrix elements of the EOM operators and of the BRST-exact operators vanish^. 

We here concentrate on diagrams (e)-(k) which are relevant for an example discussed in sect. 4: 
Diagrams (g)-(k) have already been computed in the present scheme, and the results have been 
obtained as Z ik (i, k — 1, ■ ■ ■ , n — 2) of eqs.(14)-(17) in (|TT 
renormalization constants Z ik , 
identified as 



In the singlet case, however, these 
which we denote by (Z ik ) N s here and in the following, should be 



ikJNS 



Z? k F + (-l) k nCt?Z? E G , 



(24) 

! } a c/^ a 7^- 1 ^, 
. On the other 



because T Uy E of (|18|) contains the quark-gluon-quark term i n ~ 1 S{G at11 D^ 2 ■ ■ ■ D^ 71 
which can be reexpressed by R n , k with the coefficient in the second term of (^_ 
hand, Z im , ZiE and Z mm given by eq.(18) of [|Il[] should be identified as Zf^.Zf^ and Z^. 
From the computation of diagrams (e) and (f), we obtain 

,N f 4 



X, 



FG 
IE 



n 2 CI 



n-l ' 



(25) 
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where Nj is the number of quark flavors, and I — 1, • • • , (n— 1)/2. For / = (n— l)/2 + l, ■ • ■ , n — 2, 
the replacement I — > I' — n— I — 1 should be understood. Q2~3] ) and fl2"o| ) determine Zf k F for general 
n, and the result agrees with that of ||. Similar computation including other diagrams (a)-(d) 
gives other renormalization constants. These results will be published elsewhere. 

The RG equation for the twist-3 singlet operators can be obtained in a standard way from 
As is well known, the RG equation decouples for each n. In the leading-logarithmic 



approximation (LLA), this equation is solved to give 



E«^V)» 

k=l 

n-2 



FF 



LT 



FF 



+ E« T ^V)» 



i=2 



726) 
(27) 



GF 



Ik 



Lf 



GG 

2£ 

if 



GF 2 

lm + E«^V)» 

2N f ) /3 for N c color. 



HiV, 



((RnAQ 2 ))) : 

((K, m (g 2 ))) = 
((^,(g 2 ))) = 

*!=i 

with i = 1, • • • , n — 2, I = 2, • • • , L = a s (fi 2 )/a s (Q 2 ), and 6 
Here we introduced the reduced matrix element ((0(fi 2 ))) by 

(PS-| 0^1-^-1^ | p5 ^ = _ 2 ((0( / u 2 )))<S.4 pw...p^-i] - traces, (29) 

with O = R n ,i,Rn,mi and T n ^. fl2"6|)-(|2"8D exhibit the complicated mixing characteristic of the 
higher-twist operators. On the other hand, the EOM operators R U} e, T n> E as well as the BRST- 
exact operator T n> £ decouple from the result because their matrix elements vanish ||. 

Before ending this section, we give explicit formula for the Q 2 -evolution of g2{x,Q 2 ) in the 
LLA: In the LLA, we substitute, in (§, /j 2 = Q 2 and H 2i (|, 1, a s (Q 2 )) = (e 2 )J (f - l) +0(a s ) 

for i = u,u,d,d,- ■ •, with e the charge matrix of the quarks. Therefore, the Q 2 -evolution of 
g 2 {x,Q 2 ) is governed directly by the scale dependence of the local composite operators corre- 
sponding to the moment of <f) 2 (£,Q 2 ). As shown in (^), the twist-2 contributions contained in 
92{x,Q 2 ) are determined completely by gi(x,Q 2 ) whose Q 2 -evolution is well known. On the 
other hand, for the twist-3 part g2(x,Q 2 ), we have already determined an independent basis of 
relevant operators and their scale dependence. Corresponding to the decomposition of the charge 
matrix squared by e 2 = (e 2 )l + A 3 X 3 + A 8 X 8 with A; the flavor matrices, we decompose #2(2, Q 2 ) 
into the singlet and the nonsinglet parts by g2{x,Q 2 ) = ^f(x,Q 2 ) + g2 S (x,Q 2 ). Then, from 
(D-(|13D and @, it is straightforward to see that the moment M n [g^(Q 2 )] = $ dxx n ~ l g^(x, Q 2 ) 
is given by the matrix elements of our independent operators as (see also [|IOf ) 



(e 2 ) 



J2(n-l-k)((R n , k (Q 2 )}} + 



\k=l 



n — 1 



n 



((RnAQ 2 ))) 



(30) 



The equations (^)- (|30|) completely determine the (Revolution of the singlet part. 

The moment of the flavor nonsinglet part M n [g2 S (Q 2 )] is given by the equation similar to 
fl30D with R n ,k{Q 2 ) , R n ,m(Q 2 ) and (e 2 ) replaced by the corresponding nonsinglet quantities. In 
this case, the relevant scale dependence is given by (|26|) , (|2"T|) with the replacement X FF — > (X)ns 
and ((T n j(/i 2 ))) — > 0. For the detail of the nonsinglet part, we refer the readers to e.g. JT0|, pT| . 



4 Application to the n = 3 case 

We now present an example of the Q 2 -evolution for the lowest (n = 3) moment. In this case 
there exists no three-gluon operator A • T° A , while there exists one quark-gluon-quark operator 
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A-R^ = k^G^A M fiip (see (|5 



If we neglect the contribution of the quark mass operator 



A-i?3 m , only one operator A-R% 1 contributes to the Q -evolution of the physical matrix elements 
as given by 
@: From 

eqs.(14)-(17) of 0, Zff 
where C F = (iV c 2 - 1)/2N< 
Q 2 -evolution: 



Therefore, the Q -evolution is completely determined from our results (|24|), 
he relevant renormalization constant is Z^ . We obtain, from (|2"4D , (|2"5|) and 
= (Z u )ns + 3Z™ = 1 + g 2 (C F /3 - 3C G - 2N f /3) /[8tt 2 (4 - £>)], 
Cg = iV c . This result combined with (|23|) , (p6|) and ([3~0"1) gives the 



M, 



a s (/i 2 



(3C G -C F /3+2N f /3)/b 



Several comments are in order here: (T) The result (HID coincides with that of 



(31) 



though our 



approach is quite different from those works. This fact confirms the theoretical prediction (31), 
and also demonstrates the efficiency of our method, (ii) The term 2iVy/3 would be absent from 
the exponent of ([3lD , if we considered the nonsinglet casefllC], |TTJ. Our derivation shows that this 
difference of 2iVy/3 comes from Z[^, which describes the mixing between A • R% 1 and the gluon 
EOM operator A ■ T^ E . This reveals a peculiar role played by the EOM operator: It produces 
observable effects although its physical matrix elements vanish, (iii) For N c = 3 and Nf = 4, the 
exponent of ( pT|) is 101/96, while the corresponding exponent for the nonsinglet case is 77/96. 
Thus the singlet case obeys rather stronger Q 2 -evolution compared to the nonsinglet case. Fig.^| 
shows the behavior of ( |3~1~D for the case of N c = 3, Nf — 4, and Aqcd = 0.5GeV. 



0.8 



0.6 



0.4 



0.2 




singlet 



M 3 (m(i» 



10 



100 



Q^GeV 1 ) 

Figure 2: Q 2 -evolution of g2{x,Q 2 ) for the lowest (n = 3) moment 

Conclusions 



In the present study, we have developed a manifestly covariant approach to investigate the 
flavor singlet part of g2(x,Q 2 ). We gave a thorough OPE analysis of the twist-3 flavor singlet 
operators which contribute to g2{x, Q 2 )- We derived the operator identities which relate the two- 
particle operators with the three-particle ones. We have chosen the three-particle operators as 
an independent operator's basis. To identify the renormalization constants correctly, the off-shell 
Green's functions are considered. We have shown that the EOM operators as well as the BRST- 
exact operators play an important role. As an application, we computed the Q 2 -evolution of the 
lowest (n = 3) moment. Our result confirmed the prediction obtained by different methods^, 
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The off-shell Green's functions are free from the infrared singularity coming from the collinear 
configuration, as well as from any unphysical singularities. We believe that calculating the off- 
shell Green's functions is the safest method to obtain the anomalous dimensions. 
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